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Twisted rapid passage is a type of non-adiabatic rapid passage that generates controllable 
quantum interference effects that were first observed experimentally in 2003. It is shown 
that twisted rapid passage sweeps can be used to implement a universal set of quantum 
gates Qu that operate with high-fidelity. The gate set Qu consists of the Hadamard 
and NOT gates, together with variants of the phase, 7r/8, and controUed-phase gates. 
For each gate g in Qjj, sweep parameter values are provided which simulations indicate 
will produce a unitary operation that approximates g with error probability Pe < lO"**. 
Note that all gates in Qu are implemented using a single family of control-field, and the 
error probability for each gate falls below the rough-and-ready estimate for the accuracy 
threshold Pa ~ 10'^. 

Keywords: fault-tolerant quantum computing, accuracy threshold, quantum interfer- 
ence, group-symmetrized evolution, non-adiabatic dynamics 

1 Introduction 

The accuracy threshold theorem [IHH] provides the impetus for the work presented in this 
paper. This remarkable theorem established that an arbitrary quantum computation could 
be done with an arbitrarily small error probability, in the presence of noise, and using im- 
perfect quantum gates, if the following conditions are satisfied. (1) The computational data 
is protected by a sufficiently layered concatenated quantum error correcting code. (2) Fault- 



"Corresponding author: fgaitan@lps.umd.edu 

''Permanent Address: Laboratory for Physical Sciences, 8050 Greenmead Drive, College Park, MD 20740. 



1 



2 High-fidelity universal quantum gates 



tolerant protocols for quantum computation, error correction, and measurement are used. 
(3) A universal set of unencoded quantum gates is available with the property that each gate 
in the set has an error probability Pf, that falls below a value Pa known as the accuracy 
threshold. The value of the threshold is model-dependent, though for many. Pa ^ 10"'* has 
become a rough-and-ready estimate. Thus gates are anticipated to be approaching the accu- 
racies needed for fault-tolerant quantum computing when P^ < 10^^. One of the principal 
challenges facing the field of quantum computing is determining how to implement a universal 
set of unencoded quantum gates for which all gate error probabilities satisfy P^ < 10"'*. 

In this paper numerical simulation results are presented which indicate that a class of 
non-adiabatic rapid passage sweeps, first realized experimentally in 1991 [9], and known as 
twisted rapid passage (TRP), should be capable of implementing a universal set of unencoded 
quantum gates Gu that operate non-adiabatically and with gate error probabilities satisfying 
Pe < lO^''. The gate set Gu consists of the one-qubit Hadamard and NOT gates, together 
with variants of the one-qubit vr/S and phase gates, and the two-qubit controUed-phase gate. 
The universality oi Gu was demonstrated in Ref. |T0]. This level of gate accuracy is largely 
due to controllable quantum interference effects that arise during a TRP sweep [1 H I12 ] . and 
which were observed in 2003 using NMR [13] . To find TRP sweep parameter values that yield 
such high-performance quantum gates, it proved necessary to combine numerical simulation 
of the Schrodinger dynamics with optimization algorithms that search for minima of Pe . In 
the case of the two-qubit modified controlled-phase gate, to achieve Pe < 10~*, it was also 
necessary to interleave the TRP sweep with the group-symmetrized evolution of Ref. [M] . 

The outline of this paper is as follows. Following this Introduction, Section [2] summarizes 
the essential properties of TRP, and how the numerical simulation and optimization are done. 
In Section |3] we explain how group-symmetrized evolution is incorporated into a TRP sweep, 
and then present our simulation results for each of the gates in Gu in Section |4l We close 
in Section [5] with a summary of our results, and a discussion of possible directions for future 
work. 

2 Twisted Rapid Passage 

This Section presents the essential properties of twisted rapid passage (TRP). A more detailed 
presentation of the discussion in Sections 12.11 and 12.21 appears in Refs. [10] and [12]. In an 
effort to make the present paper more self-contained, we summarize that discussion here. 
Section [2?T] introduces TRP as a generalization of adiabatic rapid passage in which the control- 
field twists in the x-y plane at the same time that its z-component undergoes a non-adiabatic 
inversion. We describe how controllable quantum interference effects arise as a consequence 
of the twisting. Section 12.21 discusses the details of the numerical simulations, and describes 
the optimization procedures used to find minima of the gate error probabilitjl? 

2.1 TRP Essentials 

To introduce TRP [HldJ], we consider a single-qubit interacting with an external control- 
field F(i) via the Zeeman interaction Hz{t) = — cr • F(i), where crj are the Pauli matrices 
{i — x,y,z). TRP is a generalization of adiabatic rapid passage (ARP) [15]. In ARP, the 
control-field F{t) is slowly inverted over a time Tq such that F(i) — atz + bic. In TRP, 
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however, the control-field is allowed to twist in the x-y plane with time-varying azimuthal 
angle 4>{t), while simultaneously undergoing inversion along the z-axis: 

F{t) ^ atz + b cos (f){t) X -I- 6 sin (/)(t) y. (1) 

Here —To/2 < t < To/2, and throughout this paper, we consider TRP with non-adiabatic 
inversion. As shown in Ref. |12j , the qubit undergoes resonance when 

h d(j) 

For polynomial twist, the phase profile <j){t) takes the form 

0„(t) = -sr. (3) 

n 

In this case, Eq. ^ has n — 1 roots, though only real- valued roots correspond to resonance. 
Ref. 11 showed that for n > 3, the qubit undergoes resonance multiple times during a 
single TRP sweep: (i) for all n > 3, when B > 0; and (ii) for odd n > 3, when B < 0. 
For the remainder of this paper we restrict ourselves to i? > 0, and to quartic twist for 
which n = 4 in Eq. ([3]). During quartic twist, the qubit passes through resonance at times 
t = 0, zty^a/hB pjj. It is thus possible to vary the time separating the resonances by varying 
the TRP sweep parameters B and a. 

Ref. [n] showed that these multiple resonances have a strong influence on the qubit 
transition probability, allowing transitions to be strongly enhanced or suppressed through a 
small variation of the sweep parameters. Ref. jl6l calculated the qubit transition amplitude 
to all orders in the non-adiabatic coupling. The result found there can be re-expressed as the 
following diagrammatic series: 

Lower (upper) lines correspond to propagation in the negative (positive) energy-level, and the 
vertical lines correspond to transitions between the two energy-levels. The calculation sums 
the probability amplitudes for all interfering alternatives |17) that allow the qubit to end up 
in the positive energy-level given that it was initially in the negative energy-level. As we have 
seen, varying the TRP sweep parameters varies the time separating the resonances. This in 
turn changes the value of each diagram in Eq. (j4]) , and thus alters the interference between the 
alternative transition pathways. It is the sensitivity of the individual alternatives/diagrams 
to the time separation of the resonances that allows TRP to manipulate this quantum inter- 
ference. Zwanziger et al. |13j observed these interference effects in the transition probability 
using NMR and found quantitative agreement between theory and experiment. It is this link 
between interfering quantum alternatives and the TRP sweep parameters that we believe 
underlies the ability of TRP to drive high-fidelity (non-adiabatic) one- and two-qubit gates. 

2.2 Simulation and Optimization Procedures 

As is well-known, the Schrodinger dynamics is driven by a Hamiltonian H{t) which causes 
a unitary transformation U to be applied to an initial quantum state lip). In this paper, it 
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is assumed that the Hamihonian H{t) contains terms that Zeeman-couple each qubit to the 
TRP control-field F{t). Assigning values to the TRP sweep parameters {a,b, B,Tq) fixes the 
control-field F{t), and in turn, the actual unitary transformation Ua applied to \ip). The task 
is to find TRP sweep parameter values that produce an applied gate Ua that approximates 
a desired target gate Ut sufficiently closely that its error probability (defined below) satisfies 
Pe < 10^**. In the following, the target gate Ut will be one of the gates in the universal set 
Gu- Since Gu contains only one- and two-qubit gates, our simulations will only involve one- 
and two-qubit systems. 

For the one-qubit simulations, the Hamiltonian Hi{t) is the Zeeman Hamiltonian Hz{t) 
introduced in Section im Ref. [T5] showed that it can be written in the following dimensionless 
form: 



Here: r = {a/b)t; X = ha/b-^; and for quartic twist, 4>4{t) — (?74/2A)t'*, with 774 — hBb^/a^. 

For the two-qubit simulations, the Hamiltonian H2(t) contains terms that Zeeman-couple 
each qubit to the TRP control-field, and an Ising interaction term that couples the two 
qubits. Alternative two-qubit interactions can easily be considered, though all simulation 
results presented below assume an Ising interaction between the qubits. To break a resonance- 
frequency degeneracy W12 ~ W34 for transitions between, respectively, the ground and first- 
excited states {El o E2) and the second- and third excited states [E^ o E4), the term 
C4\Eii{t)) {Ei{t)\ was added to H2{t). Combining all of these remarks, we arrive at the following 
(dimensionless) two-qubit Hamiltonian jlO| : 



Here: (i) bi = h'-fiBrf /2, Ui = jiBq, ji is the gyromagnetic ratio for qubit i, and i = 1,2; 
(ii) T = {a/b2)t, A = fia/b\, and 774 = fiBb^la^; and (iii) di = (wi — aj2)&2/a, c?2 = (A/a)62, 
= and ^4 — {J/a)b2, where A is a detuning parameter |10j . 

The numerical simulations assign values to the TRP sweep parameters and then integrate 
the Schrodinger equation to obtain the unitary transformation Ua produced by the resulting 
TRP sweep. Given Ua, Ut, and the initial state \tp), it is possible to work out [12] the error 
probability Pe{i') for the TRP final state |-0a) ~ Ua\4'), relative to the target final state 
IV't) = The gate error probability is defined to be the worst-case value of Pe{i')'- 



'Hi(t) = (1/A) {-TCT^ - COB (pi{T) (7 J. - SVCL (j) i{T) <7 y) . 



(5) 



^2(r) 



[-{di + d2)/2 + t/X] a\ - {d^/X) [cos04CT^ -|- sm(j)4<Ty] 
+ [-d2/2 + T/X]al -(1/A) [cos<?!)4cr^ -hsin(;()4cr^] 
-{7:d,/2)alol +C4\E4{t)){E4{t)\. 



(6) 



Pe = maxPe(V')- 



(7) 



It proves useful at this point to introduce the positive operator 




(8) 



Ref. fT2l showed that the error probability Pe satisfies the upper bound 



Pe < Tr P, 



(9) 
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where the RHS is the trace of the positive operator P. Once Ua is known, Tr P is easily 
evaluated, and so it is a convenient proxy for Pe which is harder to calculate. Tr P also has 
the virtue of being directly related to the gate fidelity 



where n is the number of qubits acted on by the gate. It is straightforward to show TO' that 



The simulations calculate Tr P, which is then used to upper bound the gate error probability 
Pe via Eq. dll). 

To find TRP sweep parameter values that yield highly accurate non-adiabatic quantum 
gates, it proved necessary to combine the numerical simulations with function minimization 
algorithms |18j that search for sweep parameter values that minimize the Tr P upper bound. 
The multi-dimensional downhill simplex method was used for the one-qubit gates, while sim- 
ulated annealing was used for the two-qubit modified controUed-phase gate. This produced 
the one-qubit gate results that will be presented in Section 14.11 However, for the modified 
controlled-phase gate, simulated annealing was only able to find sweep parameter values that 
gave Pe < 1.27 x 10"'^ [HF. To further improve the performance of this two-qubit gate, it 
proved necessary to incorporate the group-symmetrized evolution of Ref. [M] to obtain a 
modified controlled-phase gate with Pe < 10"^. In the following Section we describe how 
group-symmetrized evolution is incorporated into a TRP sweep. 

3 Group-symmetrized Evolution and TRP 

Ref. [2] introduced a unitary group-symmetrization procedure that yields an effective dy- 
namics that is invariant under the action of a finite group Q. To incorporate this group- 
symmetrization into a TRP sweep, the first step is to identify the group Q with a finite 
symmetry group of the target gate Ut, and then apply the procedure of Ref. 14 to filter out 
the ^-noninvariant part of the TRP dynamics. As the fj-noninvariant dynamics is manifestly 
bad dynamics relative to Ut, group-symmetrized TRP yields a better approximation to Ut- In 
this section we describe how group-symmetrization works, and then show how it can be incor- 
porated into a TRP sweep. The simulation results presented in Section IT2l for the two-qubit 
modified controlled-phase gate are for a group-symmetrized TRP sweep. 

3.1 Static Hamiltonian 

Consider a quantum system Q with time- independent Hamiltonian H and Hilbert space %. 
The problem is to provide Q with an effective dynamics that is invariant under a finite group 
Q, even when H itself is not t^-invariant. This symmetrized dynamics manifests as a Q- 
invariant effective propagator U that evolves the system state over a time t. Let {pi — p{gi)} 
be a unitary representation oi G onJi, and let \G\ denote the order of Q. The procedure begins 
by partitioning the time-interval (0,i) into N subintervals of duration A^^v = t/N, and then 
further partioning each subinterval into \G\ smaller intervals of duration dt^ = A.tN/\Q\. 
Let SUn = exp [— (i/?i)(5tjv-ff] denote the _ff-generated propagator for a time-interval St^, 




(10) 




(11) 
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and assume that the time to apply each pi (z G is neghgible compared to St^ (bang-bang 
limit ^19 J. In each subinterval, the fohowing sequence of transformations is apphed: 



NPi 



(12) 



The propagator U over the full time- interval (0,t) is then 



U= lim [U{Mn] 

N^oo 



N 



(13) 



Ref. ,14] showed that: 

1. for TV > 1, U{Mn) exp \-{i/n)AtNH] , where H = (l/\g\) J^fJi P^Hp^■, 

2. H is ^-invariant (viz. [H, pi] — for all pi E Q); 

3. the propagator U over (0,t) is U — exp —{i/fi)tH , which is fj-invariant due to the 
C/-invariance of H . 

The outcome of this procedure is an effective propagator U that is ^/-invariant as desired. 
3.2 Time-varying Hamiltonian 

In this subsection we show how the above procedure can be generalized to allow for a time- 
varying Hamiltonian H{t) which is the appropriate context for a TRP sweep. Although 
the generalization is straight-forward, we are not aware of any prior treatment of group- 
symmetrized evolution in the presence of a time-dependent Hamiltonian. The essential idea 
is to divide the time interval (0, t) into sufficiently small subintervals that H{t) is effectively 
constant in each. The above time-independent argument can be then be applied to each 
subinterval, yielding a ^-symmetrized propagator for that subinterval. Time-ordering the 
subinterval effective propagators then gives the effective group-symmetrized propagator for 
the entire interval (0, t). Having sketched out the basic idea, we now present the details. 

Consider a quantum system Q evolving under the action of a time-varying Hamiltonian 
H{t). We begin again by partitioning the time-interval (0,<) into JV subintervals {ti^i,ti) of 
duration Atjv = t/J\f, where i — I, . . . ,J\f. The number of subintervals is chosen sufficiently 
large that H{t) is effectively constant over each subinterval: H{t) ~ H{ti) for all t G (ti-i,ti). 
We estimate the value of Af for a TRP sweep in Section [?751 Note that for a static Hamiltonian, 
this requirement is true for all values oiAf. Since H{t) « H{ti) in the z"* subinterval (i^^i, ti), 



we can apply the symmetrization procedure of Section 13.11 to this subinterval. The result is 
the C/-invariant effective propagator 



where 



' ' i—1 



(14) 



(15) 
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is the t?-invariant effective Hamiltonian for (ti-i, U). Having the effective propagator for each 
subinterval, the effective group-symmetrized propagator U for the entire time-interval (0,i) 
is simply the time-ordered product of the U{ti,ti^i): 



U = T 



exp (^i/h dTHij)^ 



(16) 



where T indicates a time-ordered exponential, and 

1 

' ' i—l 

As in Section 13. 1[ we assume the group-symmetrizing pulses pi can be applied in a time that 
is much less than l^^tj^ /\Q\. 

3.3 TRP 

We now describe how group-symmetrized evolution can be incorporated into a TRP sweep. 
For our two-qubit simulations, the target gate is the modified controUed-phase gate 

yc,-{^^[{i' + -l)i'~{i'-<^l)^l] (18) 

which is invariant under the group Q = {I^ P , a\P , al , (j\al} . Thus \Q\ — 4, and we set 
, P4 = o\a'\. Switching over to dimensionless time, we partition the sweep 
time-interval (— to/2,to/2) into M subintervals of duration tdj^ = tq/M. We want N to 
be sufficiently large that the two-qubit Hamiltonian H2{t) is effectively constant within each 
subinterval. Since the interference effects arise from the twisting of the control field, we 
estimate the size of M by requiring that the angle swept through by the control field in a 
time Atj\f is small compared to the final twist angle 0/ = 0(to/2). Specifically, if we require 
that A(t)/(f)f < 5 X 10"'^, and noting that A(/)(t) w (j){T)Atj\f is largest at r = tq, it follows 
that 

At^f ^ A(j)/^ 

< O.OO5[</)//0(to/2)] 

< 7.5 X IQ-^. 

RecaUing that Atj^ — tq/N, and noting that tq — 120 for the two-qubit simulations (see 
Section |42|) . one finds that 

N > 1600. 

We shall see in Section HT^ that N = 2500 in the two-qubit simulations which produces an even 
smaller A0/0^, and so enhances the approximation of a constant H2{t) in each subinterval. 

The two-qubit numerical simulations for group-symmetrized TRP partition the sweep 
interval (— ro/2,To/2) into M subintervals as just described. Each subinterval (ti_i,ti) is 
further partitioned into \Q\ — 4 sub-subintervals by introducing intermediate times tij = 
ti-1 + j{Atj\f / A), with J = 1, . . . , 4. The numerical integration over the i*'' subinterval begins 
by driving the quantum state using a TRP sweep from ti^i — ti^i, at which time the group- 
symmetrizing pulse p2p\ is applied. The integration then resumes with TRP evolving the 
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state from ti^i — )■ ti^2, followed by application of P3P2 time ti^2- This alternation of TRP- 
driven evolution and group-symmetrizing pulses continues until the final time U is reached. 
The numerical integration begins at — to/2 and continues across the subintervals until the final 
time To/2 is reached. Note that the time required to apply the group-symmetrizing pulses is 
assumed to be much shorter than At_^f/4. These pulses are applied in the simulation as matrix 
operations on the state. This is consistent with our bang-bang assumption, although it pre- 
cludes a study of the effects of pulse imperfections on gate performance. This limitation will 
be removed in future work. Completion of the numerical integration yields the fj-symmetrized 
TRP propagator U corresponding to a particular choice of the sweep parameters (A,774,to) 
and system parameters (04, di, . . . , ^4). U then serves as Ua in the simulated annealing op- 
timization which returns optimized values for the sweep and system parameters, as well as 
the group-symmetrized gate Ua that best approximates the modified controlled-phase gate 
Vcp- To reduce the size of the search space, only (774,04,^1,(^4) are optimized; (tq. A, (i2, ds) 
are assigned values. The results of this simulation/symmetrization/optimization procedure 
appear in Section We shall see below that C/-symmetrized TRP yields an approximation 



to Vcp with Pe < 10-^. 



4 Simulation Results 

Here we present our simulation results for the TRP-driven approximations to the gates in 
the universal set Qu- The one-qubit gates are discussed in Section 14.11 while the two-qubit 
modified controlled-phase gate appears in Section 14.21 

4.1 One-qubit Gates 

Operator expressions for the one-qubit target gates in Qjj are: 

Hadamard : Uh = (1/a/2) ( <Tz + dx) ; 
NOT 
Modified 7r/8 
Modified phase 



Unot — 5 

K-/8 = COS {it/ 8) (Jx - sin (tt/S) ct^; 
Vp^{llV2){ax-(Ty). 



A study of the TRP-implementation of these gates was first reported in Ref. [12] . The essential 
results are included here for the reader's convenience, though space limitations do not allow 
inclusion of the unitary operator Ua associated with each gate. The interested reader can find 
them displayed in Ref. [12]. The connection between the TRP experimental and theoretical 
parameters is given in Refs. [10] and [12] for superconducting and NMR qubits, respectively. 
Table [1] presents the values for the sweep parameters A and 7^4 that produced our best 

Table 1. Simulation results for the one-qubit gates in Qu . The error probability for each gate 
satisfies P„ < Tr P. 



Gate 


A 




TrP 


T 


Hadamard 
NOT 
Modified 7r/8 
Modified phase 


5.8511 
7.3205 
6.0150 
5.9750 


2.9280 X 10-"* 
2.9277 X lO-"* 
8.1464 X lO""* 
3.8060 X lO-"* 


8.82 X lO"'' 
1.10 X 10-^5 
3.03 X 10^5 
8.20 X 10-^5 


0.9999 98 
0.9999 97 
0.9999 92 
0.9999 80 



results for TrP for each of the one-qubit gates in Qu. In all one-qubit simulations, the 
dimensionless inversion time was tq = 80.000. Since Pg < Tr P, we see that Pg < 10^"* for all 
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one-qubit gates in Qjj . Table [T] also gives the fidelity T for each gate which is obtained from 
Tr P via Eq. (fTT|) . We will discuss the robustness of these gates in Section [5l 

4-2 Modified Controlled-phase Gate 

We complete the universal set Qu by presenting our simulation results for the C/-symmetrized 
TRP implementation of the modified controlled-phase gate Vcp. In the two-qubit computa- 
tional basis (eigenstates of 



(I 

1 



\ 



\ 



-1 

1 / 



(19) 



TRP implementation of Vcp without symmetrized evolution was reported in Ref. |10| . The 
results presented there are superceded by the CJ-symmetrized TRP results presented here. 

For purposes of later discussion, note that the parameters appearing in 'H2{t) (see Eq. ©) 
fall into two sets. The first set consists of the TRP sweep parameters (A,774,to), while the 
second set (c4,c?i, . . . ,^4) consists of system parameters for degeneracy-breaking, detuning, 
and coupling. We partitioned the TRP sweep into M = 2500 pulses sequences, with each 
sequence based on the four element symmetry group Q for Vcp introduced in Section [3.3l The 
optimized parameter values A = 5.04, 774 — 3.0 x 10"'^, To = 120.00, C4 = 2.173, di — 99.3, 
d2 = 0.0, da = —0.41, and c?4 = 0.8347 produced the following two-qubit gate Ua'- 



Re {Ua) 



Im {Ua) 



( 


0.9999 98 


-0.0000 03 


-0.000015 


-0.000014 


\ 




0.0000 03 


0.9999 97 


0.0000 36 


0.0002 61 




-0.000015 


0.0000 34 


-0.9999 80 


-0.003818 




V 


-0.0000 14 


-0.0002 57 


-0.0038 38 


0.9999 81 


/ 


/ 


-0.002151 


0.0000 03 


-0.000010 


-0.0000 73 


\ 




-0.0000 03 


-0.002180 


0.000140 


-0.0003 25 






0.000010 


-0.001140 


0.001702 


0.0045 34 




V 


-0.0000 73 


-0.0003 28 


-0.0045 21 


-0.0017 78 / 



(20) 



(21) 



(i) TrP ^ 8.87 X 10"^ (u) gate fidelity Fcp = 0.9999 89; 
We see that by adding symmetrized evolution to a TRP sweep we 



Vcp, we find: 



From Ua and Ut 
and (iii) Pc < 8.87 x 10 
obtain an approximation to Vcp with Pc < 10^^. 

We see that it has been possible to use TRP sweeps to produce a high-fidelity universal 
set of quantum gates Gu, with each gate error probability falling below the rough-and-ready 
estimate for the accuracy threshold for fault-tolerant quantum computing: Pg < 10~^. 

5 Discussion 

We have presented simulation results which suggest that TRP sweeps should be capable 
of implementing the universal quantum gate set Qu non-adiabatically and with gate error 
probabilities satisfying Pc < 10"''. It is worth noting that all gates in Qjj are driven by 
a single type of control field (TRP), and that the gate error probability for all gates in 
Gu falls below the rough-and-ready estimate of the accuracy threshold Pa ^ 10~*. These 
results suggest that TRP sweeps show promise for use in a fault-tolerant scheme of quantum 
computing. 
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To achieve this high level of performance in our current formulation of TRP, some of 
the TRP parameters must be controlled to high precision. For the one-qubit gates [12], the 
critical parameter is rj^ which must be controlled to five significant figures to achieve best 
gate performance. For the modified controlled-phase gate Vcp, the critical parameters are not 
the TRP sweep parameters. Instead, for Vcp without symmetrized evolution pi);, the critical 
parameters are C4, di, and c?4 which also require five significant figure precision. However, 
when group-symmetrized evolution is added, not only is TRP able to make an approximate Vcp 
with Pg < 10^"*, but gate robustness is also improved. Specifically, di ceases to be a critical 
parameter, and C4 and d^ now only need to be controllable to four significant figures. Table [5] 
shows how Tr P varies when either C4 or c?4 is varied in the fourth significant figure, with all 

Table 2. Sensitivity of Tr P to small variation of C4 and for the two-qubit gate Vcp- AH other 
parameter values are as given in the text. 



C4 TrP 


di TrP 


2.172 6.79 X 10"^ 

2.173 8.87 X 10^5 

2.174 7.73 X 10"^ 


0.8346 1.52 X 10"^ 
0.8347 8.87 x 10"^ 
0.8348 1.52 X lO"^ 



other parameters held fixed. Thus, adding group-symmetrized evolution improves both the 
accuracy and robustness of the TRP approximation to Vcp- Note that four significant figure 
precision corresponds to 14-bit precision which can be realized with present-day arbitrary 
waveform generators (AWG) [20]. On the other hand, the current precision requirements 
for all one-qubit gates in Qij lie beyond the reach of existing commercially available AWGs. 
Unfortunately, group-symmetrized evolution cannot be used to improve the robustness of the 
one-qubit TRP gates. It is possible to show that Ut = a ■ cr, the only one-qubit unitary 
operators that commute with Ut are the identity and a multiple of Ut- Thus the only symmetry 
group available that does not include Ut is the trivial group whose sole member is the identity. 
Some other means must be found to improve the robustness of the TRP approximations to the 
one-qubit gates in Qu. Two approaches are currently under study based on: (i) the Hessian 
of our cost function Tr P; and (ii) quantum optimal control theory. 

In previous work [21] we have studied a number of forms of polynomial, as well as periodic, 
twist. To date we have found that quartic twist provides best all-around performance when it 
comes to making the gates in Qu. Although we do not at present have arguments that explain 
why this is so, ongoing work based on quantum optimal control provides a framework with 
which this question can be studied. This represents an important direction for future work. 

Finally, Refs. [10]- [12] have shown how TRP sweeps can be applied to NMR, atomic, 
and superconducting qubits. We note that TRP-generated quantum gates should also be 
applicable to spin-based qubits in quantum dots as such qubits also Zeeman-couple to a 
magnetic field. 

Acknowledgements 

We thank: (i) F. Nori, RIKEN, and CREST for making our visit to RIKEN possible; (ii) RIKEN 
for access to the RIKEN Super Combined Cluster on which the simulations incorporating 
group-symmetrized evolution were done; and (iii) T. Howell HI for continued support. 

References 



Li and Gaitan 11 



1. D. Aharanov and M. Ben-Or, Fault-tolerant computation with constant error, in Proceedings of 
the Twenty-Ninth ACM Symposium on the Theory of Computing, 176 (1997). 

2. A. Y. Kitaev, Quantum computation algorithms and error correction, Russ. Math. Surv. 52, 1191 
(1997). 

3. A. Y. Kitaev, Quantum error correction with imperfect gates, in Quantum Communication, Com- 
puting, and Measurement (Plenum Press, New York, 1997), pp. 181-188. 

4. D. Gottcsman, Stabilizer codes and quantum error correction, Ph. D. thesis, CaUfornia Institute 
of Technology, Pasadena, CA (1997). 

5. E. Knill, R. Laflamme, and W. H. Zurek, Resilient quantum computation. Science 279, 342 (1998). 

6. E. Knill, R. Laflamme, and W. H. Zurck, Resilient quantum computation: error models and thresh- 
olds, Proc. R. Soc. Lond. A 454, 365 (1998). 

7. J. Preskill, Reliable quantum computers, Proc. R. Soc. Lond. A 454, 385 (1998). 

8. F. Gaitan, Quantum error correction and fault-tolerant quantum computing (CRC Press, Boca 
Raton, FL 2008). 

9. J. W. Zwanzigcr, S. P. Ruckcr, and G. C. Chingas, Measuring the geometric component of the 
transition probability in a two-level system, Phys. Rev. A 43, 3232 (1991). 

10. R. Li, M. Hoover, and F. Gaitan, High-fidelity universal set of quantum gates using non-adiabatic 
rapid passage, Quantum Info. Comp. 9 290 (2009). 

11. F. Gaitan, Temporal interferometry: a mechanism for controlling qubit transitions during twisted 
rapid passage with a possible application to quantum computing, Phys. Rev. A 68 052314 (2003). 

12. R. Li, M. Hoover, and F. Gaitan, High-fidelity single-qubit gates using non-adiabatic rapid passage. 
Quantum Info. Comp. 7 594 (2007). 

13. J. W. Zwanziger, U. Werner- Zwanziger, and F. Gaitan, Non-adiabatic rapid passage, Chem. Phys. 
Lett. 375 429 (2003). 

14. P. Zanardi, Symmetrizing evolutions, Phys. Lett. A 258 77 (1999). 

15. A. Abragam, Principles of nuclear magnetism (Oxford University Press, New York 1961). 

16. F. Gaitan, Berry's phase in the presence of a non-adiabatic environment with an application to 
magnetic resonance, J. Mag. Reson. 139 152(1999). 

17. R. P. Feynman and A. R. Hibbs, Quantum mechanics and path integrals, (McGraw-Hill, New 
York, 1965). 

18. W.H. Press et al., Numerical Recipes (Cambridge University Press, New York 1992). 

19. L. Viola and S. Lloyd, Dynamic suppression of decoherence in two-state quantum systems, Phys. 
Rev. A 58, 2733 (1998). 

20. For example, the Tektronix AWG5000B arbitrary waveform generator provides 14-bit vertical 
resolution. 

21. R. Li and F. Gaitan, Controlling qubit transitions through quantum interference during non- 
adiabatic rapid passage, Optics and Spectroscopy 99, 257 (2005). 



